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摘 X. 积 雪 是 冰冻 圈 中 较为 活跃 的 因子 ,对 气候 环境 变化 敏感 ,其 变化 影响 着 全 球 气候 和 水 文 的 


变化 。 积 雪 和 覆盖 日 数 (SCD)、 降 雪 


ui 


始 时 间 (SCOD ) 和 融雪 
量 平衡 的 主要 因素 。 使 用 MODIS 无 云 积 雪 产 品 提取 了 叶 尔 羌 河流 域 2002 年 7 月 一 2018 年 6 月 还 


始 时 间 (SCMD) 是 影响 地 表 物 质 和 能 


日 积 雪 和 覆盖 率 (SCP), 基 于 像 元 计算 了 SCD、SCOD 和 SCMD, 系统 地 分 析 了 其 空间 分 布 与 变化 特征 ， 
并 探讨 了 其 变化 的 原因 及 积 雪 面 积 的 异常 变化 与 ENSO 的 联系 。 结 果 表 明 :(1) 研究 时 段 内 ,流域 


UTE E 


积 呈 微弱 减少 趋势 ,与 气温 呈 显 著 负 相关 ,与 降水 呈 显 著 正 相关 ;2002 一 2018 年 ,SCP 


随 海拔 的 升 高 呈 明 显 的 线性 增加 趋势 ( 尼 =0.92 .P<0.01) ); 各 海拔 高 度 带 最 大 SCP 出 现 的 月 份 大 臻 
随 海 拔 的 上 升 往 后 推迟 ,最 小 SCP 出 现 月 份 无 显著 变化 (集中 在 8 月 ) ,海拔 4000 m 以 下 ,春季 的 
SCP 小 于 冬季 ,海拔 4000 m 以 上 ,春季 的 SCP 大 于 冬季 。(2) SCD.SCOD 和 SCMD 有 明显 的 海拔 梯 
度 ,在 流域 内 ,从 东北 至 西南 ,呈现 出 SCD 增加 ,SCOD 提前 ,SCMD 推迟 的 特征 ;变化 趋势 上 ,流域 
91.9% 的 区 域 SCD 表现 为 减少 ,65.6% 的 区 域 SCOD 有 往 后 推迟 的 趋势 ,77.4% 的 区 域 SCMD 表现 出 


提前 的 趋势 。(3) 2006、.2008 年 和 2017 年 积 雪 履 盖 面 积 异 常 偏 大 ,而 在 2010 年 则 蜡 常人 篇 小 ,其 原因 


可 能 是 ENSO 影响 了 积 雪 的 变化 (4) 以 喀 喇 昆仑 为 主 的 高 海拔 地 区 ,包括 帕 米 尔 高 原 东 部 的 部 分 


地 区 ,其 SCD .SCOD 和 SCMD 分 别 表 现 出 增加 、 提 前 


走低 以 及 降水 量 的 增加 有 关 。 


和 推迟 的 趋势 ,这 种 变化 与 其 春秋 温度 的 持续 


关键 词 : 叶 尔 羌 河流 域 ; 积 雪 覆盖 日 数 ;， 降雪 开始 时 间 ; 融雪 开始 时 间 
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积 雪 作 为 冰冻 圈 中 最 活跃 的 自然 元 素 之 一 ,以 
其 特殊 的 物理 特征 (如 高 反照 率 、 高 融化 潜 热 和 低 
导热 率 ) 和 广泛 的 地 理 分 布 影响 着 地 气 系统 的 能 量 
交换 ,水 循环 和 气候 变化 一 。 积 雪 面 积 的 减少 会 引 
起 地 表 反 照 率 的 降低 ,进而 影响 下 垫 面 吸收 太阳 短 
波 辐射 ,导致 地 表 温 度 和 地 表 热 通 量 的 变化 3。 积 
雪 对 气候 变化 的 高 度 敏 感性 ,使 其 在 气候 系统 变 暖 
的 背景 下 减少 显著 。IPCC 的 最 新 报告 指出 ,几乎 所 
有 地 区 ,特别 是 在 低 海 拔 地 区 , 积 雪 持续 时 间 平 均 
10 a 减少 5 d, 男 外 ,包括 积 雪 在 内 的 冰冻 圈 内 自 
然 元 素 的 减少 已 经 改变 了 大 多 数 相 关 自 然 灾害 发 
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生 的 频率 .强度 和 位 置 5。 

近年 来 ,遥感 技术 以 其 多 尺度 、 多 时 相 、 多 谱 
段 、 多 层次 等 特点 为 相关 研究 者 开展 高 海拔 山区 的 
积 雪 研究 提供 了 相对 更 为 优质 的 数据 来 源 呈 ,已 有 
不 少 学 者 对 不 同 地 区 的 积 雪 时 空 变化 及 其 对 气候 
因子 的 响应 进行 了 研究。 其 中 一 些 研 究 基 于 多 源 
遥感 数据 在 大 尺度 上 对 积 雪 进行 了 时 空 变 化 分 析 ， 
发 现 北 半球 的 积 雪 面积 自 1922 年 以 来 都 呈现 出 减 
少 的 趋势 ; 亚 欧 大 陆 春季 融雪 开始 时 间 在 1972 一 
2006 年 显著 提前 ;中 国 积 雪 徐 盖 日数 在 2000 一 
2014 年 ,夏季 呈 减 少 趋势 ,而 在 冬季 、 春 季 和 秋季 都 
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呈现 增加 的 趋势 " ;青藏 高 原 中 部 积 雪 持续 日 数 在 
2001 一 2014 年 呈 增 加 趋势 "…。 也 有 一 些 研究 者 在 
流域 尺度 上 对 积 雪 变 化 进行 了 研究 ,如 Zheng SE 
对 新 疆 玛 纳 斯 河流 域 积 雪 的 垂直 分 布 特征 及 其 与 
气温 的 关系 进行 了 人 研究 ;Ahmad 等 5 对 印度 河 支流 
奇特 拉 尔 河流 域 积 雪 面 积 的 时 空 变 化 趋势 进行 了 分 
析 , 发 现 奇特 拉 尔 河流 域 的 积 雪 范 围 在 整个 流域 尺 
度 以 及 各 个 高 度 带 都 旦 显著 增加 趋势 ;Hasson SE 
All Tahir 等 都 研究 了 印度 河 子 流 域 阿 斯 托 河流 
域 和 洪 扎 河流 域 的 积 雪 变 化 , Hasson 等 观察 到 阿 斯 
托 河流 域 和 洪 扎 河流 域 的 积 雪 在 减少 ,而 Tahir 等 的 
研究 表明 阿 斯 托 河 流域 和 洪 扎 河流 域 的 积 雪 面 积 
略 有 增加 。 影 响 积 雪 变 化 的 气象 因子 中 ,气温 和 降 
水 是 影响 积 雪 年 内 分 配 与 年 际 变化 以 及 空间 分 布 
的 两 大 主要 气象 因子 ,已 有 不 少 研究 者 就 积 雪 
对 气候 变化 的 响应 进行 了 研究 ,Xiong 等 "通过 对 
高 亚洲 地 区 近 35 a 融雪 开始 时 间 序 列 的 分 析 及 其 与 
平均 气温 的 比较 , 发现 平均 气温 与 融雪 开始 日 期 具 
有 较 强 的 相关 性 , 除 喀 喇 昆仑 山 和 昆仑 山西 部 部 分 
地 区 外 ,1988 一 2015 年 高 亚洲 地 区 的 融雪 开始 时 间 
逐渐 提前 ; 曾 小 自 等 所 利用 MODIS 积 雪 数 据 对 土 库 
曼 斯 坦 山区 的 积 雪 进行 了 监测 ,发 现 积 雪 面 积 随 气 
温 的 降低 、 降 雨量 的 减少 而 增加 ; 姜 康 等 中 研究 发 
现 ,温度 是 影响 积 雪 日 数 、 积 雪 初 日 和 积 雪 终日 等 
积 雪 特 性 参数 变化 的 重要 因素 。 

叶 尔 羌 河流 域 属于 新 疆 干旱 区 ,流域 内 的 喀 喇 
昆仑 山 和 东 帕 米尔 高 原 是 我 国 现代 冰川 作用 的 主 
要 区 域 之 一 ,尤其 喀 喇 昆仑 山 是 亚洲 山岳 冰川 最 集 
中 的 地 区 2 ,冰雪 融 水 是 叶 尔 羌 河流 域 重 要 的 淡水 
来 源 ”。 在 全 球 变 暖 背景 下 ,北半球 的 积 雪 积 累 期 
普遍 缩短 ,消融 期 普遍 延长 , 积 雪 覆盖 面积 普遍 减 
少 , 但 是 各 区 域 由 于 气候 和 地 形 等 条 件 的 不 同 , 积 
雪 的 积累 和 消融 过 程 也 有 所 不 同 ,研究 叶 尔 羌 河 流 
域 的 降雪 开始 时 间 融雪 开始 时 间 和 积 雪 覆盖 日 数 
的 时 空 变化 特征 对 进一步 理解 其 积 雪 的 积累 和 消 
融 过 程 具有 重要 意义 。 因 此 ,本 文采 用 MODIS 无 云 
积 雪 产品 ,研究 了 叶 尔 羌 河流 域 积 雪 帮 盖 范围 的 变 
化 以 及 降雪 开始 时 间 、 融 雪 开 始 时 间 和 积 雪 覆盖 日 
数 的 时 空 分 布 和 变化 特征 ,并 探讨 了 其 变化 的 原因 
及 积 雪 面积 异常 变化 与 ENSO 的 联系 ,以 期 为 叶 尔 
羌 河 流域 的 洪水 泥石流 和 雪灾 预警 等 研究 提供 新 
依据 ,为 冰雪 水 资源 的 合理 利用 提供 新 参考 。 


1 研究 区 概况 


人 研究 区 位 于 青藏 高 原 西 北部 边缘 的 叶 尔 羌 河 
FLER, TPF 35°28’ ~39°45'N FU 73°41 ~78°47'E Z M] 
(图 1)。 流 域 总 面积 98241 kn ,南部 为 喀 喇 昆仑 山 
和 西 昆 仑 山 ,西北 部 为 帕 米 尔 高 原 , 东 依 塔 克 拉 玛 
于 沙漠 边缘 ,东北 部 为 出 前 冲积 平原 ,地 势 西 南 高 
东北 低 , 平 均 海 拔 约 3339 m。 流 域内 主要 的 山脉 为 
喀 喇 昆仑 山 ( 北 坡 ), 按 特定 的 地 理 位 置 和 地 形 条 件 
大 致 分 为 山区 与 平原 两 大 气候 区 ,山区 为 高 原 高 寒 
干旱 半 干 旱 气候 ,平原 属 干 旱 半 干旱 气候 ,全 年 干 
旱 少 十 六 。 流 域内 发 育 着 丰富 的 冰川 (图 1 中 的 冰 
川 数 据 采 用 中 国 第 二 次 冰川 编目 数据 汪 ) ,冰川 融 
水 占 河流 多 年 平均 径流 量 的 64% ,降雨 和 积 雪 融 水 
补给 占 13.4%, 地 下 水 补给 占 22.6% ,地 面 观测 数 
据 显 示 叶 尔 羌 河 流域 多 年 平均 气温 在 3.6~12.7 °C, 
年 平均 降水 量 在 57.3~78.8 mm”。 受 西风 气流 控 
制 ,其 降水 主要 集中 在 冬 春 两 季 , 流域 高 大 的 山 势 
和 高 山 冷 储 作 用 使 得 降水 主要 以 固态 水 的 形式 降 
落下 来 ,高 大 山体 截留 高 空 水 汽 的 能 力 又 使 得 高 
上 的 降水 随 海拔 逐渐 递增 ,为 积 雪 的 持续 发 育 提供 
了 较为 丰富 的 物质 来 源 ”。 


À 


图 例 


C] Landsat OLI 边 界 


高 程 /m 
- 8513 


= 0 80km 
1089 I= 


图 1 研究 区 概况 图 
Fig. 1 Sketch map of the study area 


2 数据 与 方法 


2.1 数据 

2.1.1 积 雪 数 据 本 文 的 积 雪 数据 源 于 青藏 高 原 地 
X MODIS 逐日 无 云 积 雪 产品 光 ,空间 分 辨 率 为 500 
m, 时 间 范 围 是 2002 年 7 月 一 2018 年 6 月 ,由 中 国 科 
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学 数据 中 心 提供 (http://www.csdata.org/en/p/)。 与 地 
面 观测 数据 相 比 ,卫星 获取 的 观测 数据 具有 空间 覆 
R .时间 分 辩 率 高 等 优点 ,尤其 在 人 口 稀少 的 由 
米尔 高 原 和 喀 喇 昆仑 山地 区 ,观测 站 点 十 分 稀少 ， 
卫星 积 雪 产 品 为 研究 未 观测 地 区 的 积 雪 特征 提供 
了 重要 的 数据 集 。 

为 验证 MODIS 积 雪 数据 在 研究 区 的 精度 ,利用 
空间 分 辨 率 为 30m ÉI Landsat OLI 数 据 [ 来 源 于 美国 
地 质 勘 探 局 (USGS) , https://landsatlook.usgs.gov/ 
viewer.html] 通 过 snowmap 算 法 提取 积 雪 面积 ,并 使 
用 决定 系数 (R) .偏差 (Bias) , 均 方 根 误差 (RMSE )、 
相对 均 方 根 误差 (CRMSE) 来 评价 MODIS 积 雪 数据 
在 研究 区 的 表达 情况 。 验 证 数据 选取 的 时 段 为 
2014 年 7 月 1 日 一 2015 年 6 月 30 日 整个 积 雪 年 , 包 
括 图 1 所 示 的 研究 区 中 所 有 云 覆 盖 率 <10% 的 影像 
(行列 号 为 150033、149033、149034、148034 和 
148035) ,共计 32 景 。 经 统计 ,MODIS 积 雪 数 据 与 
landsat OLI 提 取 的 积 雪 数据 一 致 性 较 好 ,能 较 好 地 
反映 研究 区 的 积 雪 时 空 分 布 特征 (图 2)。 
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TER’ Bias RMSE FI rRMSE 分 别 表示 决定 系数 、 俩 差 . 均 方 根 
误差 和 相对 均 方 根 误差 。 
图 2 MODIS 积 雪 产品 的 精度 验证 
Fig.2 Validation of MODIS snow cover product 


2.1.2 和 气象 数据 由 于 研究 区 域内 气象 站 点 非常 有 
限 ,为 研究 气温 和 降水 的 空间 差异 对 积 雪 分 布 及 变 
化 的 影响 ,本 文采 用 格 网 气象 数据 ,以 避免 站 点 数 
据 代 表 的 空间 范围 有 限 带 来 的 影响 。 数 据 来 源 于 
Climatic Research Unit (CRU) (https://crudata.uea.ac. 
uk/cru/data/hrg/cru_ts_4.03/) ,采用 TS 数据 集 , 版 本 为 
v4.03 ,为 逐 月 气象 数据 ,空间 分 辨 率 为 0.5?x0.5"” 。 


此 外 ,本 文 以 2002 年 7 月 一 2018 年 6 月 流域 内 莎 车 
和 喀什 气象 站 点 观测 的 气温 和 降水 数据 [中 国 地 面 
气候 资料 日 值 数据 集 (V3.0) ,来 源 于 国家 气象 信息 
FF» Chttp///data.ema.en/) ] 验 证 同期 CRU 的 气温 和 
降水 数据 ,气温 数据 之 间 的 相关 系数 高 达 0.99 ,降水 
数据 的 相关 数据 在 0.63~0.67 之 间 ,说 明 CRU 和 气温 和 
降水 数据 在 研究 区 的 精度 较 高 ,能 满足 研究 需要 。 
2.1.3 DEM &4& DEM 数 据 采用 SRTM (Shuttle ra- 
dar topography mission) ,数据 版 本 为 V.003 ,来 源 于 
美国 地 质 勘 探 局 (USGS) (https://Ipdaac.usgs.gov/ 
products/srtmgl 1v003/) ,空间 分 辨 率 为 30 m。 本 文 利 
JH DEM 数据 研究 不 同 海拔 高 度 的 积 雪 特 征 ,为 与 积 
雪 数 据 相 匹配 ,将 其 重 采样 和 重 投影 与 积 雪 数据 一 
样 的 分 辨 率 和 投影 。 

2.1.4 ENSO 指数 MEI 表征 厄尔尼诺 -南方 涛 动 
(ENSO) 的 指数 MEI(Multivariate ENSO Index) (数据 
版 本 为 V2) 由 NOAA 提供 (https://www.esrl.noaa.gov/ 
psd/enso/mei/)。 相 比 于 其 他 的 ENSO 指数 ,如 SOI 和 
Nimo 3.4 SST, MEI 更 加 完整 有 效 地 描述 了 ENSO 
现象 ,能 较 好 的 反映 海 气 看 合 的 性 质 , 且 不 易 受到 
数据 更 新 的 偶然 误差 造成 的 影响 中 。 本 文 利 用 
2002—2018 年 的 MEI 指 数 探究 积 雪 异常 与 大 气 环 
流 异常 的 联系 。 

2.2 方法 

2.2.1 积 雪 参数 的 提取 FAG mK (Snow cover 
percentage, SCP) 代 表 相 对 积 雪 和 覆盖 面积 ,是 人 研究 区 
积 雪 和 窗 盖 面积 占 人 研究 区 总 面积 的 百分比 , 积 雪 窗 盖 
日 数 (Snow covered days, SCD ) 指 一 个 积 雪 年 中 有 积 
oy eT AY KB, BS FF UR EY TH] (Snow cover onset 
dates, SCOD ) 代 表 一 个 积 雪 年 中 首次 降雪 的 日 期 ， 
融雪 开始 时 间 (Snow cover melting dates, SCMD ) 则 
表示 一 个 积 雪 年 中 积 雪 开始 融化 的 日 期 ,本 研究 基 
于 Wang 等 ”提出 的 算法 计算 每 个 积 雪 年 的 降雪 开 
始 时 间 、 融 雪 开 始 时 间 。 为 便于 积 雪 特 征 的 研究 ， 
选取 当年 7 月 1 日 一 次 年 6 月 30 日 为 一 个 积 雪 年 。 
本 文 利 用 逐日 积 雪 禾 盖 产品 计算 了 叶 尔 羌 河流 域 
每 日 的 SCP, 并 基于 像 元 ,计算 每 个 像 元 中 的 SCD、 
SCOD 以 及 SCMD ,具体 计算 公式 如 下 : 


P, 
SCP (1) 


n 


SCD = V (S) (2) 


i=l 


SCOD = Fdl - SCDI (3) 
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SCMD = Fd2 +SCD2 (4) 
式 中 : P 代表 研究 区 中 积 雪 有 覆盖 的 面积 ; P 代表 人 研 
究 区 的 总 面积 ;n 代 表 一 个 积 雪 年 内 包含 的 天 数 ; S, 
为 1 表示 为 积 雪 日 ,为 0 表示 为 非 积 雪 日 ;Fd1 为 冷 
季 中 国定 的 某 一 天 ,鉴于 叶 尔 羌 河 流域 的 降雪 一 般 
多 发 于 9 一 12 月 ,Fdl 取 12 月 1 日 , 即 Fd1=335 d zX 
336 d;SCD1 为 每 个 积 雪 年 中 9 月 1 日 一 12 月 1 日 的 
积 雪 窗 盖 日 数 ;Fd2 为 冷 季 中 国定 的 某 一 天 ,鉴于 流 
域 的 融雪 一 般 多 发 于 2 一 5 月 ,Fd2 取 2 月 1 日 , 即 
Fd2=32 d,SCD2 为 每 个 积 雪 年 中 2 月 1 日 一 5 月 1 日 
的 积 雪 覆盖 日 数 。 
2.2.2 时 序 分 解 ” 由 于 本 文 所 用 的 积 雪 数据 和 气象 
数据 均 为 时 序数 据 , 要 明确 时 序数 据 的 具体 变化 特 
征 , 就 需要 分 解 时 序 以 获得 精细 信息 。STL(Season- 
al-trend decomposition procedure based on Loess ) 作 为 
一 种 常用 的 时 序 分 解 方法 ,广泛 应 用 于 气象 和 水 文 
相关 的 研究 ,其 将 时 间 序 列 Y 分 为 趋势 分 量 7 周期 
分 量 (季节 分 量 ) S 及 余 项 R ,以 精细 化 地 表征 时 序 
数据 在 不 同 层面 的 变化 特点 下 。 对 于 时 间 序 列 
Y ,其 i 时 刻 的 SIL 分解 可 表示 为 式 (5)"。 
Y,=7,+S,+R, (=1, 2, +++, N) (5) 
研究 过 程 中 主要 使 用 STL BET RR A 
温 、 降 水 等 时 序数 据 的 分 解 , 以 分 析 其 各 分 量 的 特 
征 ,发 气 有 效 信息 。 
2.2.3 趋势 检验 方法 本文 利 用 Mann-Kendall 非 参 
数 检验 方法 和 Theil-Sen Median #4344) Br TE f oU E. 
HE EXT SCD ,SCOD 和 SCMD HEFT EAA) Hr 和 趋势 
Aim, 


3 结果 与 分 析 


3.1 积 雪 范 围 时 空 变化 特征 

图 3 展示 了 叶 尔 羌 河流 域 2002 一 2018 年 多 年 
平均 SCP 的 年 内 变化 情况 ,由 图 可 知 积 雪 从 9 月 开 
始 积累 ,缓慢 增长 至 11 月 ,从 12 月 开始 快速 增长 ， 


7 8 9 10 1112 12 34 5 6 
Att 


图 3 2002 一 2018 年 叶 尔 羌 河 流域 积 雪 和 覆盖 率 
(SCP) 的 年 内 变化 
Fig. 3 Annual cycle of SCP in the Yarkant River 
Basin during 2002—2018 


间 序 列 进 行 分 解 , 获 得 其 各 自 的 趋势 项 .季节 项 及 
残 差 项 (图 4)。 从 SCP 的 季节 项 可 以 发 现 , 积 雪 覆 
盖 面 积 的 季节 变化 特征 十 分 明显 ,2、3 月 为 SCP 的 
峰值 期 ,7、8 月 为 谷 值 期 。 通 过 对 SCP 的 时 间 序 列 
做 线性 趋势 拟 合 , 得 出 其 变化 的 斜率 为 -0.0164, 说 
HH 2002—2018 年 流域 的 积 雪 和 覆盖 面积 呈 微 弱 减 少 
趋势 。SCP 趋 势 项 代表 积 雪 禾 盖 面积 的 总 体 变 化 趋 
势 , 从 叶 尔 羌 河流 域 2002 一 2018 年 SCP 的 变化 趋势 
可 以 发 现 ,从 2002 年 至 2006 年 , 积 雪 面 积 呈 现 出 上 
升 的 趋势 ,2006 年 至 2007 年 , 随 着 气温 的 上 升 和 降 
水 的 减少 , 积 雪 面积 又 减 ,然后 至 2012 年 , 积 雪 面 积 
呈 微 弱 的 波动 上 升 趋势 ,至 2013 年 , 随 着 气温 的 上 
升 , 积 雪 面 积 再 度 观 减 ,而 后 缓慢 上 升 至 2017 年 ,到 
2018 年 , 随 着 气温 的 上 升 和 降水 的 减少 , 积 雪 面积 
再 次 又 减 , 可 以 发 现 SCP 与 降水 呈 总 体 一 致 的 变化 
趋势 ,而 与 气温 呈 相 反 的 变化 趋势 。 进 一 步 将 流域 
每 月 平均 SCP 与 每 月 平均 气温 和 降水 做 相关 性 分 
析 ,发 现 SCP 与 降水 呈 显 著 正 相关 ,相关 系数 为 
0.61 ,与 气温 呈 显 著 负 相关 ,相关 系数 为 -0.73( 表 
1) ,说 明 流 域 积 雪 禾 盖 面 积 受气 温和 降水 的 影响 ， 
且 对 气温 的 敏感 性 大 于 对 降水 的 敏感 性 。SCP 的 残 


直至 次 年 2 月 达到 峰值 (35.94% ) ;之 后 , 随 着 春季 气 
温 的 回暖 , 积 雪 开 始 加 速 消 融 , 至 8 月 SCP 降 至 最 低 
值 ; 此 外 ,SCP 在 冬季 的 波动 较 大 ,其 中 12.1 月 和 2 
月 的 波动 尤其 大 ,夏季 SCP 的 波动 相对 较 小 ,其 中 7 
月 和 8 月 的 波动 最 小 。 

为 研究 积 雪 覆盖 面积 的 年 际 变化 特征 ,采用 
2.2.2 部 分 的 方法 对 流域 SCP 以 及 气温 和 降水 的 时 


差 项 可 以 反映 流域 积 雪 有 覆盖 面积 波动 的 特征 ,是 极 
端 降雪 事件 最 有 力 的 证 据 , 在 一 定 程 度 上 可 以 反映 
极端 气候 事件 ,由 图 4 中 SCP 残 差 项 的 变化 趋势 可 
Jj, , 2006 ,2008 4F Fil 2017 4E AY SCP FFE fa. , Mi 2010 
年 的 SCP 则 异常 偏 小 ,这 可 能 与 这 些 年 份 的 气候 异 
BAX. 

由 图 Sa 可 知 ,SCP 表 现 出 明显 的 海拔 梯度 , 随 
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图 4 2002 一 2018 年 叶 尔 羌 河 流域 积 雪 覆盖 率 (SCP) 以 及 气温 、 降 水 变化 
Fig. 4 Variations of SCP, temperature, and precipitation in the Yarkant River Basin during 2002—2018 


表 1 叶 尔 羌 河流 域 积 雪 覆 盖 率 (SCP) 与 气温 、 
降水 的 相关 系数 
Tab.1 Pearson correlation coefficients between SCP 
and temperature, precipitation in the Yarkant River 
Basin on monthly basis 
气温 降水 
TECH sic -0.73* 0.61" 
注 :** 表 示 在 0.01 水 平 上 显著 相关 。 


海拔 的 升 高 呈 明 显 的 线性 增加 趋势 (R=0.92、P< 
0.01) , 积 雪 和 覆盖 面积 以 5200~5250 m 海 拔高 度 带 为 
中 心 向 低 海拔 和 高 海拔 呈正 态 分 布 (图 5b)。 图 6 为 
不 同 海拔 高 度 带 多 年 平均 SCP 的 年 内 变化 ,可 以 发 
现 海拔 越 高 ,SCP 的 年 内 波动 越 小 ;海拔 大 于 6500 m 
的 区 域 各 月 的 SCP 都 接近 100% ;海拔 小 于 4000 m 
的 高 度 带 春季 的 SCP 小 于 冬季 ,而 大 于 4000 m 的 高 
度 带 春季 的 SCP 大 于 冬季 ,这 可 能 是 因为 在 海拔 较 
高 的 区 域 冬 季 积 雪 升 华 、 风 吹雪 导致 积 雪 发 生 了 再 
分 布 , 从 而 造成 积 雪 大 量 损失 。 此 外 , 受 西风 影 
响 ,海拔 4000~5000 m 和 5000-6000 m 的 SCP TE A 
春 两 季 较 高 。 各 海拔 高 度 带 最 大 SCP 出 现 的 月 份 大 
致 随 着 海拔 的 上 升 往 后 推迟 ,海拔 2000~3000 m 以 
及 3000~4000 mm 高 度 带 最 大 SCP 出 现在 2 月 , 随 着 海 
拔 的 上 升 ,海拔 4000~5000 m 以 及 5000~6000 m ij 
度 带 最 大 SCP 出现 的 月 份 依次 往 后 推迟 为 3 月 和 5 
月 ;与 此 同时 , 除 海拔 2000~3000 m 以 及 3000~4000 
m 高度 带 的 SCP 分 别 在 5 月 和 7 月 达 最 小 值 ,其 余 海 
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图 5 叶 尔 羌 河流 域 积 雪 履 盖 率 (SCP)(a) 、 积 雪 和 覆盖 面积 和 
流域 面积 随 海拔 变化 (b) 
Fig. 5 Variations of the annual mean SCP (a), snow cover 
area and basin area under different altitude in the Yarkant 
River Basin (b) 


IEH 4000 m 以 上 的 高 度 带 最 小 SCP 均 出 现在 8 月 。 
3.2 积 雪 参 数 时 空 变化 特征 
由 叶 尔 羌 河流 域 2002 一 2018 年 多 年 平均 积 雪 
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Fig.6 Mean monthly variations of SCP for different 


elevation zones in the Yarkant River Basin during 2002—2018 


日 数 的 空间 分 布 可 以 看 出 (图 7a) , 受 地 势 和 气候 等 
因素 的 影响 ,SCD 大 致 呈现 出 由 东北 向 西南 递增 的 
特征 ,海拔 越 高 ,SCD 越 高 (图 8) ,例如 , 喀 喇 昆仑 、 
西 昆仑 以 及 帕 米 尔 高 原 的 高 海拔 冰川 区 ,其 SCD 大 
于 240 d, 而 小 于 60 d 的 SCD 低 值 区 主要 分 布 在 流域 
中 海拔 较 低 的 河谷 和 山 前 冲积 平原 。SCOD 的 分 布 
随地 势 的 变化 则 表现 出 由 东北 向 西南 逐渐 提早 的 
变化 特点 ,高 海拔 冰川 区 的 SCOD 普遍 较 早 ,9 月 初 
就 开始 出 现 积 雪 ,主要 分 布 在 喀 喇 昆仑 山 以 及 帕 米 


尔 高 原 东部 和 西 昆仑 山 的 部 分 地 区 ; 随 着 海拔 的 降 
低 ,SCOD 逐渐 推迟 ,至 流域 中 海拔 较 低 的 河谷 和 下 
游 地 区 ,11 月底 至 12 月 初 才 出 现 积 雪 ,在 空间 分 布 
上 ,SCOD 有 了 明显 的 海拔 梯度 ,海拔 越 高 ,SCOD 越 早 
(图 8)。 而 SCMD 则 表现 出 由 东北 向 西南 逐渐 往 后 
推迟 的 特征 ,SCMD 的 这 种 分 布 特征 主要 是 受 流域 
地 势 特征 的 影响 , 积 雪 消 融 最 早 在 流域 中 海拔 较 低 
的 河谷 和 下 游 地 区 ,主要 开始 于 2 月 上 旬 , 随 着 海拔 
的 上 升 ,融雪 初 日 逐渐 后 推 , 至 高 海拔 的 冰川 区 , 融 
雪 一 般 开 始 于 5 月 上 旬 和 中 旬 , 同 SCD 和 SCOD 一 
样 ,SCMD 也 有 比较 明显 的 高 度 梯 度 , 随 海拔 的 升 高 
SCMD 逐渐 往 后 推迟 。 

变化 趋势 上 (图 7b~c) , 2002—2018 年 HEAR IE 
河流 域 91.9% 的 区 域 SCD 表现 为 减少 的 趋势 ,其 中 
显著 减少 的 区 域 主 要 分 布 在 由 米尔 和 西 昆仑 山 的 
部 分 地 区 ,增加 的 区 域 主要 位 于 喀 喇 昆仑 山 和 帕 米 
尔 高 原 的 少 部 分 地 区 ,这 与 孙 燕 华 等 中 观察 到 的 帕 
米尔 高 原 的 积 雪 变化 一 致 。 从 SCOD 的 变化 趋势 来 
看 (图 7e~f) ,流域 大 部 分 地 区 的 SCOD 略 有 推迟 ,其 
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7 2002 一 2018 年 叶 尔 羌 河流 域 积 雪 履 盖 日 数 (SCD ) 降雪 开始 时 间 (SCOD ) Fil 


融雪 开始 时 间 (SCMD ) 的 空间 分 布 .变化 趋势 及 显著 性 


Fig.7 Annual mean distribution, variation trend, and its significance of SCD, SCOD, 
and SCMD in the Yarkant River Basin during 2002—2018 
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图 8 2002 一 2018 年 叶 尔 羌 河 流域 积 雪 覆盖 日 数 (SCD)\ 降 

雪 开始 时 间 (SCOD) 和 融雪 开始 时 间 (SCMD ) 随 海拔 变化 

Fig. 8 Variations of SCD, SCOD, and SCMD with elevation 
in the Yarkant River Basin during 2002—2018 


rp 65.589685 K E SCOD 有 往 后 推迟 的 趋势 ,但 仅 2% 
的 区 域 SCOD 显著 推迟 ,主要 分 布 在 西 昆 仑 和 帕 米 
尔 高 原 的 东部 ,提早 的 区 域 主要 位 于 喀 喇 昆仑 山 和 
帕 米 尔 东部 的 部 分 地 区 ,尤其 是 喀 喇 昆仑 山 的 部 分 
地 区 表现 出 显著 提前 的 趋势 ,提前 速率 >2 deat. Tit 
域 中 77.4% 的 区 域 SCMD 表现 出 提前 的 趋势 (图 
7h~i) ,但 其 中 仅 2.06% 的 区 域 SCMD 显著 提前 ,主要 
分 布 在 西 昆 仑 山 以 及 帕 米 尔 东部 ,以 喀 喇 昆仑 为 主 
的 高 海拔 地 区 ,包括 帕 米尔 高 原 东部 的 部 分 地 区 的 
融雪 开始 时 间 呈 现 出 推迟 的 趋势 。 


4 讨论 


4.1 大 气 环流 与 积 雪 变化 

只 雪 变 化 主要 是 大 尺度 大 气 环 流 或 气候 强迫 
作用 的 结果 ” ,例如 西风 环流 .季风 ,北极 涛 动 和 
厄尔尼诺 -南方 涛 动 等 。 叶 尔 羌 河流 域 位 于 青藏 高 
原 西 北部 ,其 积 雪 的 分 布 特征 受 西风 环流 运动 路 径 
以 及 水 汽 输送 方向 的 影响 。 西 风 带 高 空冷 槽 通过 
青藏 高 原 西 北部 上 空 时 , 受 地 形 阻 流 作用 ,出 现 分 
文 , 北 文 绕 过 天 山 , 南 支 沿 青 藏 高 原 西 侧 南下 , 叶 尔 
芜 河 流域 西部 的 由 米尔 高 原 以 及 西南 部 的 喀 喇 昆 
仑 山 正 处 于 南 广西 风 的 上 升 运动 区 ,因此 ,降水 较 
多 ,在 高 海拔 地 区 降水 以 积 雪 的 形式 降落 ,使 得 其 
积 雪 覆盖 日 数 较 大 ; 当 气 流 仆 坡 翻越 流域 西部 和 西 
南部 的 山区 向 东北 前 行 时 ,所 剩 水 汽 较 少 , 使 得 流 
域 山 区 的 北部 虽 海 拔 较 高 ,但 因 离 水 汽 较 远 ,降水 
较 少 ,因此 , 积 雪 履 盖 日 数 也 较 少 “”。 

积 雪 异 常 与 大 气 环 流 异常 密切 相关 , 当 环 流 形 
势 出 现 异 常 变化 时 ,就 会 直接 影响 天 气 和 气候 ,其 


至 引发 极端 气候 事件 ,造成 降雪 和 积 雪 异常 。 
有 关 研 究 表 明 ,大 气 环 流 遥 相关 通过 改变 温度 、 降 
水 以 及 大 气 的 动力 和 热力 作用 ,从 而 引起 积 雪 的 变 
化 ”” 引 , 其 中 ,以 厄尔尼诺 -南方 涛 动 (ENSO) 琐 相关 
最 为 典型 。ENSO 是 热带 太平 洋 海 气 相互 作用 的 结 
果 , 具 有 2~7 a 的 准 周期 ,对 气候 的 季节 -年 际 变化 
起 着 重要 的 作用 ”“。 从 4.1.1 部 分 的 分 析 中 可 知 ， 
2006 、2008 年 和 2017 年 叶 尔 羌 河流 域 积 雪 覆盖 面积 
异常 偏 大 ,而 2010 年 的 积 雪 覆盖 面积 则 异常 偏 小 。 
为 探究 其 异常 的 原因 ,使 用 表征 ENSO 的 指数 MET 
和 流域 SCP 的 残 差 项 以 及 气温 和 降水 的 残 差 项 进行 
分 析 。 从 图 9a~p 可 以 看 出 ,2006.2008 年 和 2017 年 
的 MEI 值 连续 数 月 小 于 零 ,与 此 同时 ,这 些 年 份 冬 春 
季 的 气温 偏 小 ,降水 偏 大 ,其 中 尤 以 2008 年 最 为 典 
型 ,说 明 这 些 年 份 可 能 发 生 了 ENSO 冷 事件 , 即 拉 尼 
娜 ;而 2010 年 冬季 的 MEI 值 连续 数 月 大 于 零 ,上 且 气 
温 偏 大 ,降水 偏 小 ,说 明 2010 年 可 能 发 生 了 ENSO BE 
事件 , 即 厄 尔 尼 诺 。 因 此 ,推测 叶 尔 羌 河流 域 积 雪 覆 
盖 面 积 的 异常 增 大 ( 减 小 ) 可 能 是 由 于 拉尼娜 ( 厄 尔 
尼 诺 ) 事 件 通 过 海 气 作用 的 遥 相 关 , 对 流域 的 气候 产 
生 了 一 定 影响 ,从 而 引起 了 降雪 和 积 雪 的 异常 。 
4.2 气候 变化 与 积 雪 变 化 

积 雪 变化 与 气候 变化 密切 相关 ,气温 和 降水 作 
为 主要 的 气象 因子 ,研究 其 空间 分 异 和 变化 特征 及 
与 积 雪 变 化 的 联系 对 揭示 积 雪 分 布 和 变化 的 成 因 
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图 9 2002 一 2018 年 叶 尔 羌 河 流域 SCP 残 差 项 与 MEI 时 间 
序列 (a) 以 及 气温 与 降水 残 差 项 时 间 序列 (b) 
Fig. 9 Variations of MEI, residual component of SCP (a), 


temperature and precipitation (b) in the Yarkant River 
Basin during 2002—2018 
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及 与 气候 的 相互 关系 据 有 重要 意义 。 已 有 不 少 研 
究 发 现 相对 于 降水 ,气温 是 影响 积 雪 变化 的 主要 因 
素 ”“"。 本 文通 过 对 比 积 雪 覆 盖 率 与 气温 和 降水 的 
相关 系数 ,发 现 这 一 地 区 积 雪 覆盖 率 与 气温 的 相关 
性 大 于 与 降水 的 相关 性 ,表明 气温 是 影响 叶 尔 疙 河 
流域 积 雪 面积 变化 较为 主要 的 因素 。 受 地 势 的 影 
响 , 叶 尔 疙 河流 域 的 年 平均 气温 表现 出 由 西南 向 东 
北 递增 的 空间 变化 特征 ,年 均 降 水 量 的 分 布 则 大 致 
呈现 出 由 西南 向 东北 递减 的 特征 ,以 喀 喇 昆仑 ( 北 
坡 )、 西 昆仑 和 帕 米 尔 高 原 东 部 为 主 的 高 海拔 山区 


人 研究 发 现 , 叶 尔 羌 河流 域 的 SCD 与 其 春季 降水 总 量 
呈 显 著 正 相关 , 与 秋季 降水 总 量 呈 微弱 正 相 关 ( 表 
1) ,说 明 春 季 和 秋季 降水 量 的 增加 会 促进 积 雪 的 积 
累 , 延 长 SCD。 因 为 叶 尔 羌 河 流域 的 降雪 主要 开始 
于 秋季 , 而 融雪 主要 发 生 在 气温 接近 于 冰点 的 春 
季 , 所 以 秋季 的 降水 量变 化 和 平均 气温 很 大 程度 上 
决定 着 SCOD 的 变化 ,而 春季 降水 量 的 变化 和 春季 
平均 气温 则 决定 着 SCMD 的 变化 。 通 过 将 SCOD 和 
SCMD 与 秋季 、 春 季 降 水 总 量 和 秋季 ,春季 平均 气温 
做 相关 分 析 , 发 现 SCOD 与 秋季 降水 总 量 呈 负 相 关 ， 


的 年 平均 气温 在 零度 以 下 ,年 均 降水 量 总 量 则 在 
200 mm 以 上 (图 10a~p)。 因 此 ,气温 和 降水 的 这 种 
空间 分 异 特征 使 得 叶 尔 羌 河 流域 的 积 雪 和 覆盖 日 数 
也 表现 出 由 西南 向 东北 递减 的 特征 。 

正如 Peng 等 中 所 指出 的 ,即使 气温 升 高 , 当 气 
温 低 于 冰点 时 ,降水 加 剧 也 会 导致 积 雪 深度 增加 。 


40°N F (2 
39°N 
38°N 
37°N 


年 平均 气温 /*C 


36°N 


35°N 1 1 1 1 1 L 
T4°E 76*E 78°E 


40°N + (c) 
39°N 
38°N 
37°N 


秋季 平均 气温 /*C 


36°N 


35°N 
74°E 76°E 78°E 


40°N L (e) 
39°N 
38°N 
37°N 


春季 平均 气温 /"C 


36?N 


35?N 1 | 1 1 1 1 
T4°E 76*E T8°E 


与 秋季 平均 气温 呈正 相关 ,SCMD 与 春季 降水 总 量 
呈正 相关 ,与 春季 平均 气温 呈 负 相关 , 即 秋季 降水 
量 的 增加 和 气温 的 减少 会 促使 SCOD 提前 ,春季 降 
水 总 量 的 增加 和 和 气温 的 减少 会 促使 SCMD 往 后 推迟 
( 表 2)。 又 如 图 10c-d 所 示 ,2002 一 2018 年 位 于 叶 
尔 羌 河流 域 西南 部 的 喀 喇 昆仑 山 和 西北 部 的 帕 米 
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图 10 2002 一 2018 年 叶 尔 羌 河 流域 多 年 平均 气温 (a) 和 多 年 平均 降水 总 量 (b) 空 间 分 布 秋季 平均 气温 (c) 和 
秋季 降水 总 量变 化 趋势 (d) 以 及 春季 平均 气温 (e) 和 春季 降水 总 量 (f) 变 化 趋 热 


Fig. 10 Spatial distributions of annual mean temperature (a) and precipitation (b), variation trends of average temperature (c) 


and total precipitation in autumn (d), and that of average temperature (e) and total precipitation in spring (f) in the 
Yarkant River Basin during 2002—2018 
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表 2 叶 尔 芜 河流 域 积 雪 才 盖 日 数 (SCD ) 降雪 开始 时 间 (SCOD ) 和 融雪 开始 时 间 (SCMD ) 与 气温 \ 降 水 的 相关 系数 


Tab.2 Pearson correlation coefficients of SCD, SCOD and SCMD with temperature, precipitation in the 


Yarkant River Basin on monthly basis 


秋季 降水 总 量 秋季 平均 气温 春季 降水 总 量 春季 平均 气温 
SCD 0.24 -0.26 0.43" -0.20 
SCOD -0.40 -0.35 -0.07 
SCMD 0.02 0.15 -04T 


* 


注 :* 表 示 在 0.05 水 平 上 显著 相关 。 


尔 高 原 东 部 秋季 降水 总 量 呈 现 出 微弱 的 增加 趋势 ， 
在 其 秋季 平均 气温 小 于 零 的 背景 下 ,降水 的 增加 促 
进 了 其 积 雪 深度 的 增加 ,从 而 使 得 其 SCOD 表现 出 
提前 的 趋势 ;同样 ,如 图 10e~f 所 示 , 流 域 中 喀 喇 昆 


仑 和 帕 米 尔 东部 春季 降水 总 量 表现 为 增加 , 且 春 季 


平均 气温 小 于 零 , 低 温和 降水 量 的 增加 使 其 积 雪 次 
度 增 加 ,从 而 促使 其 SCMD 往 后 推迟 。 


5 结论 


(1) 2002—2018 年 ,流域 积 雪 覆盖 面积 呈现 出 
微弱 的 减少 趋势 , 旦 与 气温 和 降水 的 变化 关系 显 
著 , 与 气温 呈 负 相关 ,与 降水 呈正 相关 ;年 内 变化 
上 , 积 雪 覆盖 面积 在 8 月 最 小 ,2 月 最 大 ;空间 变化 
上 ,SCP 随 海拔 的 升 高 呈 明 显 的 线性 增加 趋势 (及 = 
0.92 ,P«0.01) , 积 雪 有 覆盖 面积 则 以 5200~5250 m FIR 
高 度 带 为 中 心 向 低 海拔 和 高 海拔 呈正 态 分 布 ,海拔 
4000 m 以 下 ,春季 的 积 雪 面积 小 于 冬季 ,海拔 4000 
m 以 上 ,春季 的 积 雪 面 积 大 于 冬季 ;此 外 ,各 海拔 高 
度 带 最 大 SCP 出 现 的 月 份 随 着 海拔 的 上 升 往 后 推 
iR ,最 小 SCP 出 现 月 份 无 显著 变化 (集中 在 8 月 )。 

(2) SCD( 随 海拔 的 升 高 增加 )、SCOD( 随 海拔 
的 升 高 往 前 提早 ) 和 SCMD (海拔 的 升 高 往 后 推迟 ) 
有 了 明显 的 海拔 梯度 ;在 流域 内 ,从 东北 至 西南 ,呈现 
出 SCD 增 加 ,SCOD 提前 ,SCMD 推迟 的 特征 ;研究 时 
段 内 ,流域 91.9% 的 区 域 SCD 表现 为 减少 ,65.58% 的 
区 域 SCOD 有 往 后 推迟 的 趋势 ,77.4% 的 区 域 SCMD 
表现 为 提前 。 

(3) 2006 、2008 年 和 2017 年 积 雪 覆盖 面积 异常 
偏 大 ,而 在 2010 年 则 异常 偏 小 ,其 原因 可 能 是 ENSO 
事件 引起 了 降雪 和 积 雪 的 异常 。 

(4) 以 喀 喇 昆仑 为 主 的 高 海拔 地 区 ,包括 帕 米 
尔 高 原 东 部 的 部 分 地 区 ,其 SCD、SCOD 和 SCMD 分 
别 表现 出 增加 、 提 前 和 推迟 的 趋势 ,这 种 变化 与 其 


春秋 温度 的 持续 走低 以 及 降水 量 的 增加 有 关 。 
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Spatiotemporal variation of snow cover in the Yarkant River 
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YI Ying", LIU Shiyin", ZHU Yu", WU Kunpeng" 
(1. Yunnan Key Laboratory of International Rivers and Transboundary Eco-security, Kunming 650500, Yunnan, China; 


2. Institute of International Rivers and Eco-Security Yunnan University, Kunming 650500, Yunnan, China) 


Abstract: Snow cover is one of the active elements in the cryosphere sensitive to climate changes. Its change can 
affect the global climate and hydrology changes, and its surface mass and energy balance are affected by the 
snow-covered days (SCD), snow cover onset dates (SCOD), and snow cover melting dates (SCMD). In this study, 
we employed the MODIS daily cloud- free snow cover product to calculate variations of snow cover percentage 
(SCP), SCD, SCOD, and SCMD from July 2002 to June 2018 in the Yarkant River Basin, south Xinjiang, China 
and then analyzed their distribution and spatiotemporal variations. We discussed the cause of the variation and the 
relation between the abnormal change of snow cover and ENSO. We found the following results: (1) In the study 
period, the snow cover area exhibited a slight decrease trend and had a significant negative correlation with 
temperature and a positive correlation with precipitation. SCP showed a clear linear increase trend with 
increasing altitude (R? —0.92, P«0.01). When the maximum SCP appears in each altitude zone, the month is 
generally delayed with the rise in altitude. In contrast, when the minimum SCP appears, the month has no 
obvious change with altitude (concentrated in August). SCP in spring is less than that in winter when below the 
4000 m altitude. On the other hand, SCP in spring is greater than that in winter when above the 4000 m altitude. 
(2) Clear elevation gradients were found for SCD, SCOD, and SCMD. The spatial distributions of SCD, SCOD, 
and SCMD are gradually increased, delayed, and advanced, respectively, in a northeast-to-southwest direction in 
the basin. In the basin, from 2002 to 2018, SCD decreased to 91.9%, SCOD delayed to 65.6%, and SCMD 
advanced to 77.4%. (3) In 2006, 2008, and 2017, the snow cover area was abnormally large, whereas, in 2010, it 
was abnormally small, which may be affected by ENSO. (4) High mountains concentrated along the Karakoram, 
including Parts of the eastern Pamir Plateau, experienced a lengthened SCD, an advanced SCOD, and a delayed 
SCMD related to the continuous low temperature and the increase of precipitation in spring and autumn. 


Key words: Yarkant River Basin; snow covered days; snow cover onset dates; snow cover melting dates 


